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(54) Process and system for management of electrical power supply 



(57) A computer implemented tool forecasts the 
spot price of electric power in a deregulated market and 
the amounts of power that may be traded in this market. 
The computer implemented tool makes these forecasts 



at different delivery points, providing probabilistic distri- 
butions for spot prices and trading, to allow proper risk 
management of power supply within the electricity net- 
work. 
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Description 

The present invention generally relates to a computer process and system for use in an electricity network to 
manage power supply by forecasting spot-market prices of electric power and trading transactions at different delivery 

s points based on various technical and commercial criteria. 

Electricity is an essential part of our lives. If we took a moment to think of the services that would not exist if 
electricity were not available, we would be amazed. Almost every aspect of our modern lives involves electric power, 
from light bulbs and television sets to hospitals and automobile industries. Although we are used to having power 
whenever we need it, the processes and systems involved in delivering electricity require careful planning and sophis- 

10 ticated mathematical and forecasting models, in order to allow for various technical, environmental, and commercial 
considerations. Examples of technical considerations include the operating characteristic of power stations (which may 
not be able to respond instantly to demand changes), and capacitity limitations of the transmission network; an example 
of an environmental consideration is the weather, which has a strong impact on consumer demand for electric power; 
and an example of a commercial consideration is the action of a competitor (eg the amount and price of electricity 

15 which it is prepared to supply). 

Lately, due to rising costs of energy and discrepancies in its price among different regions of the country, the legal 
framework surrounding the electric-power industry in the USA has changed (a process known as deregulation). This 
change has opened the door for utilities to compete with each other and against independent suppliers regardless of 
their geographic location. Although this change will benefit the consumer, utilities are going to face a highly unpredict- 

20 able market and will need to make tough decisions regarding power generation and delivery. The current picture of a 
single utility controlling the market in a specific region will soon disappear. Instead, there will be power producers who 
sell their product to a power pool; and power suppliers who will buy power from this pool and in turn sell it to their 
customers. Although the full picture of the power industry after deregulation is not yet known, it is clear that utilities 
need to prepare themselves for an open market in which buying and selling power are to be considered when scheduling 

25 their generating units. The main reason behind such deregulation is to reduce the high price of electric energy. Initial 
steps towards deregulation were taken in 1978 with the passage in the USA of the Public Utilities Regulatory Policy 
Act. This act encouraged non-utility generation and required utilities to buy power from independent generators. The 
Energy Policy Act of 1992 took deregulation a step further by mandating open access to the transmission system for 
wholesalers. 

30 Currently, electricity is sold as a service that is delivered to specified points. For example, one typically expects to 

receive electric power via a meter outside the house. We pay for this service regardless of its producer or which power 
lines it followed. That is, an electricity bill indicates the total usage of electricity in kilowatt hours (KWH) and the service 
price per KWH without incorporating any other details into the pricing scheme. Deregulation in the USA is changing 
this picture by unbundling the electric power into generation and transmission. In the future, one will pay a production 

35 cost and a transmission fee. There will be several power suppliers from whom electric power may be purchased. 
Suppliers may have different pricing mechanisms. For instance, there might be a discount for using power off-peak 
periods or for signing a long-term contract with the supplier. Power producers will compete with each other to minimize 
their costs so that they can sell their product to more customers and maximize their profit. 

The hope is that deregulation will result in cheaper prices and play a part in improving the economy by encouraging 

40 investment in electric utilities. The size of the electric industry is expected to grow after deregulation as was the case 
with the telecommunications industry. The telecommunications industry's revenue shot up from $81 billion to $170 
billion within ten years of deregulation. 

A pitfall of deregulation is that the load on a utility system is becoming increasingly unpredictable. The reason is 
that trading transactions can change the load pattern significantly. For example, some utilities may sell more than 30% 

45 of their power generation to other utilities on certain days. Demand and supply in the market are functions of volatile 
electricity prices which in turn depend on highly unpredictable elements such as regional weather conditions and fuel 
prices. 

There are two components for any competitive environment: physical and commercial. In the case of electricity, 
the physical component already exists in the form of transmission lines and all consumers have physical access to the 

50 market. The commercial component consists of the ability of buyers and suppliers to execute trades over the physical 
network. One should also ensure that anyone who is interested in joining this market has nondiscriminatory access to 
the physical network. In this environment, suppliers and consumers have the freedom to contract with each other under 
negotiated terms and conditions. To ensure the delivery of the power traded, they need a mechanism to ensure the 
physical delivery of the electricity. Also, they may need to rely on a spot market to buy or sell additional electricity. 

55 Although the new environment in the USA encourages competition, power transmission will remain regulated there 

for the time being in order to maintain a reliable system. The transmission lines in each state or region will be controlled 
by an independent entity called an Independent System Operator or ISO. One of the ISO responsibilities is to settle 
financially with the parties involved in transmitting electric power The transmission cost depends on the proximity of 
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the supplier and the congestion of the transmission lines as well as other operational factors. To maintain a reliable 
system, an ISO announces, twenty-tour hours in advance, the load forecast on the system and asks interested suppliers 
to submit bids; i.e., amount of power to be delivered at each point and pricing information. The load forecast is given 
for each delivery point within the physical system. Bids are submitted to the Power Exchange or PX which then holds 
an auction to determine which suppliers to buy power from. For each delivery point, the PX chooses the least expensive 
source of power, then the second least expensive, and so on. The process is repeated until the forecasted load on the 
system is met. Note that such a process ensures that the customers' demand for electricity is satisfied at a minimal 
cost. The bidding cost of last supplier chosen in the bidding process is the spot-market price for power at this delivery 
point. The PX performs this process for each delivery point in the system and at each hour of the day 

After agreeing on the production schedule, the ISO studies the schedule to make sure that it is feasible; i.e., it 
does not exceed the capacity of the transmission lines. Transmission charges are a function of the congestion of the 
power lines. These charges are calculated by ISO and billed to the different parties involved in trades. Note that sup- 
pliers who submit a bid with a high price may end up not selling any of their production. On the other hand, selling 
power at a low price may not create enough revenue for a generator. Given that no one knows in advance the amount 
of power that competitors may bid for, the electric-power market will become more uncertain and risky. The hope is 
that competition will force utilities to minimize their production costs which will, in turn, benefit consumers and the 
economy in general. This model also encourages individuals to invest money in the electric-power market. For instance, 
if the charges across a certain power line are relatively high due to congestion, an investor may build an alternative 
power line in the hope of benefitting from the high-traffic in that region. The result, of course, is a reduction in the 
transmission cost in that region and an increase in the reliability of the system. 

It will be appreciated that the USA is not alone in introducing deregulation. For example, in the UK, the formerly 
public electricity utility (a monopoly) has been split at privatisation into competing generators, and multiple distributors. 
Again, whilst such deregulation is expected to lead to lower prices, it significantly adds to the complexity of operating 
in the electricity network. In particular, it will be recognised that in the face of deregulation, the proper management of 
electric power supply requires accurate forecasting in order to avoid technical and/or commercial failure (the former 
represented by the failure of the network as a whole to meet consumer demand, the latter by the economic demise of 
a power generator or distributor). 

Accordingly, the invention provides a computer process for use in an electricity network to manage power supply 
by forecasting spot-market prices of electric power and trading transactions at different delivery points, said method 
comprising the steps of: 

inputting maximum generating capacities of each utility competing within the market, price functions at which a 
utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost and transmission 
capacity in different power lines of the electrical system; 

using probabilistic distributions of electric load and the price function of each utility and transmission cost and 
transmission capacity between different utilities to build probabilistic distributions for the spot-market prices and 
for electric trades in the system; and 

outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning horizon 
and for a plurality of delivery points to assist in said supply management. 

Such a forecasting process provides the technical and commercial information necessary for producing or distrib- 
uting power in order to achieve a satisfactory operation of the electricity network, and the provision of output probabilistic 
distributions in particular helps in effective risk managment. 

It is preferred that the process further comprises the step of defining different scenarios wherein the maximum 
generating capacity of each utility changes under said different scenarios, and also further comprises the step of using 
the input weather forecast to create a probabilistic distribution of load forecasts for each utility of interest. 

The input price functions may be given implicitly, for example by approximating a price function of a utility as a 
quadratic function which is convex, positive and increasing against certain data, or they may be given explicitly as a 
cost function of the generating units. 

In a preferred embodiment, the probabilistic distribution for the spot-market prices and for electric trades in the 
system are created by sampling from electric load, price functions and cost distributions to create a sequence of de- 
terministic models, with the cost of operating an electric utility system being minimised under each sampled set of price 
functions, transmission costs and electric loads. 

Further in the preferred embodiment, the probabilistic distributions for spot prices and trades are determined by a 
bucketing process comprising the steps of: 

creating ranges of spot prices; 

counting the number of spot prices that fall within each range and normalizing the counted numbers to obtain a 
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probability for each range; 

and generating a histogram ot the probabilities of each range as an approximation of the probability distribution 
of the marginal cost at each delivery point at a given time. 

s The invention further provides a system for use in an electricity network to manage power supply by forecasting 

spot-market prices of electric power and trading transactions at different delivery points, said system comprising: 

means for inputting maximum generating capacities of each utility competing within the market, price functions at 
which a utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost and 
10 transmission capacity in different power lines of the electrical system; 

means for using probabilistic distributions of electric load and the price function of each utility and transmission 
cost and transmission capacity between different utilities to build probabilistic distributions for the spot-market 
prices and for electric trades in the system; and 

means for outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning 
is horizon and for a plurality of delivery points to assist in said supply management. 

Viewed from another aspect the invention provides a computer power trading and forecasting process for fore- 
casting spot-market prices of electric power and trading transactions at different delivery points comprising the steps of: 

20 inputting maximum generating capacities of each utility competing within the market, price functions at which a 

utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost transmission 
capacity in different power lines of the electrical system; 

using probabilistic distributions of electric load and the price function of each utility and transmission cost and 
transmission capacity between different utilities to build probabilistic distributions for the spot-market prices and 
25 for electric trades in the system; and 

outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning horizon 
and for a plurality of delivery points to enable a power broker to effectively manage risk in buying and selling electric 
power with different suppliers and consumers. 

30 in a preferred embodiment, the probabilistic distributions of the transmission costs and capacities of power lines 

are created from historical data, and the probabilistic distributions for the spot -market prices and for electric trades in 
the system are created by sampling from electric load, price functions and cost distributions to create a sequence of 
deterministic models. Preferably, the process further comprises the steps of minimizing the cost of operating an electric 
utility system under each sampled set of price functions, transmission costs and electric loads; and storing a resulting 

35 output. The sampling from electric load, transmission costs and cost distributions is performed using a random number 
generator. 

Thus the computer implemented process described herein forecasts the spot-market prices and the trading trans- 
actions at different delivery points, providing the decision maker with probabilistic distributions for spot prices and 
trading to allow efficient risk management. The user of the process can be (1 ) a utility or an independent power supplier 

40 who is interested in using the trading and spot market information to decide which generating units to commit, how to 
price power, and the level of generation needed; or (2) a power broker who is interested in buying and selling power 
(but not producing) in order to maximize revenue. The computer implemented process uses the weather forecast in 
addition to other market information and provides probabilistic distributions for the power traded and spot-market prices, 
rather than a single estimate for the spot price and the power traded. An example of the output of the computer imple- 

45 mented process is shown in Table 1 below; the information therein can be used for hedging against uncertainty, thereby 
allowing risk to be acceptably contained. 



Table 1 
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12.00 


-50 
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Table 1 (continued) 



$ 


MWH 


Pr. 


13.00 


-30 


5 


13.00 


-20 


10 


17.50 


40 


5 


17.50 


50 


5 


18.00 


50 


15 


18.00 


60 


10 


18.00 


70 


5 



is The first column in the above table represents the spot market price. The second column contains the amount of power 
that may be traded. Note that for the same spot price, one may buy or sell different amounts of electric power. The 
likelihood of a transaction is given in the last column. 

It will be appreciated that it is desirable to use all available historical, competitive, and weather information to 
produce the most accurate output, although the data that is actually available for different utilities may differ from one 
region to another. For example, some utilities are required, by the power pool, to post the status of their generating 
units by the end of the week. These requirements differ significantly between power pools. Also, data that are currently 
public, such as hourly loads, may not be available in the near future. Accordingly, the tool can be operated in a flexible 
manner in accordance with the amount of information that is in fact available. 

A preferred embodiment of the invention will now be described in detail by way of example only with reference to 
the following drawings: 



35 



45 



55 



Figures 1 A and 1 B are block diagrams illustrating the relation of the power trading and forecasting tool to its inputs 
and outputs; 

Figure 2 is a flow diagram illustrating the process of approximating the price function of a utility by a quadratic 
function; 

Figures 3A and 3B are graphs demonstrating how the cost function of a generator can be approximated using a 
linear function; 

Figure 4 is a graph showing an approximate cost function of a utility constructed by using the linear approximations 
of Figure 3; 

Figure 5 is a flow diagram of the simulation used in the process; 

Figure 6 is a flow diagram of a process that generates a random sample; and 

Figure 7 is a flow diagram showing the bucketing process used to create a probabilistic distribution using the output 
of simulation runs. 

This application is related in subject matter to U.S. Patent Application Serial No. 08/869,561 (IBM Docket 
YO997-032) filed June 5, 1997, by Samer Takriti for "RISK MANAGEMENT SYSTEM FOR ELECTRIC UTILITIES". 
This earlier application is incorporated herein by reference, and a copy of this earlier application is available on the file 
of the present application. 

Referring now to the drawings, and more particularly to Figures 1 A and 1 B, there are shown in block diagram form 
the relationship of the power trading and forecasting tool and its inputs and outputs. The tool can have different sets 
of mput parameters, and two are described here by way of example. The first example is shown in Figure 1 A. In this 
case, we are given the historical data of selling prices and amounts of power sold by each utility. The tool 11 receives 
inputs from a plurality of utilities 12 t to ^2 n . These inputs include for each utility prices at which it sold power in the 
past and the corresponding amounts of power sold. In addition, the tool 11 receives as an input the transmission costs 
and the maximum transmission capacity for different lines in the transmission network 1 3. The output of tool 11 is the 
probabilistic distribution for trades and spot prices 14. 

The second set of input parameters is illustrated in Figure 1 B. They include the operating cost curves of the different 
generators 12-, to 12 n , estimates of the marginal profits, and electric load forecasts. The output is similar to that of 
Figure 1 A. The inputs and outputs of the tool 11 are described in more detail below. 

System's Inputs 

Since the details of the deregulated market may vary from country to country, and so cannot be completely defined, 
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the input parameters, data available about utilities and their generators, are not necessarily known in advance. To 
avoid this difficulty, we describe the input parameters in general terms so that our tool can handle any electric power 
market. The input to our system consists of (1) the generating capacity and price functions, ff t , of utilities competing 
in the market, (2) the load forecast in the regions of interest, and (3) the transmission costs, cj^ fJ charged by the ISO 
5 or other equivalent body to transmit electricity from node / to node j at time period r and the transmission capacity, 
-j , • Note that all of our notations contain the superscript s. This is the scenario index, a concept that is described 
in 'further detail in the above referenced US patent application and which shall be clarified through the following dis- 
cussion. Here is a detailed description of the input parameters. 

1. The maximum generating capacity, Gf v of each utility, /, competing within the market. Note that this capacity 
to is not necessarily the total generating capacity of all generators owned by a utility. For example, some units may be 

unavailable due to maintenance during certain days. Also, depending on the season and availability of fuel, the capacity 

may change. The scenario index, s, allows a certain utility to have different possible capacities. The capacity could 

also change depending on rover the planning horizon, usually one week or 168 periods. For example, the capacity 

over the weekends is lower than regular days. 
*5 2. The price functions, ff t at which utility / is willing to sell its power at time t under sbenario s The value 

returned by ff t is assumed to include the operating cost plus the profit margin of a utility. Note that ff t may be given 

(a) implicitly or (b) explicitly. 



The Implicit Case 



20 



We may not know the exact cost structure of a utility but may be able to approximate its operating cost by a 
quadratic function. For instance, we may observe that a utility, /, priced its power at time f and under a certain load at 
Xf f $/MWH. By observing the price offered by a utility for its power over time, one may construct an approximation of 
the cost structure of this utility. The simplest way to do so is to collect the different data points of price versus load and 
25 fit a quadratic curve using least squares. Our experience indicates that the quadratic approximation is a good one for 
most utilities. Of course, if one observes a certain profile for the cost structure of a utility, one may consider a different 
function to be fit. 

The problem of "guessing" the price function of a utility is not a simple one in a competitive environment. For 
example, one may observe a cost of $15 per MWH for a certain utility under a 12,000 MWH load on one day while 

30 observing an operating cost of $18 per MWH for the same utility under a similar load in another day. The reason for 
this discrepancy might be a generating unit that is going through maintenance, a failure of a generator in that utility, a 
large error in the load forecast which resulted in shortage of generation, or a shortage of certain fuel which may have 
forced the utility to use a more expensive fuel. Regardless of the reason, building the price function f? t of a utility /may 
be impossible if one has to use a single value for the cost under a given load. To handle this difficulty, our model permits 

35 the use of more than one price function, f? f , at any time period t. For a given time period f, we allow the price function 
to be any of the functions fj r -ff t where s/is the number of different functions available at time t. That is, we model the 
cost function of a utility / at each time period t by a set of cost functions. The index s is read as scenario. One hopes 
that using a set of functions reflects accurately the uncertainty an observer may have regarding the operating cost of 
a utility. For example, one scenario may represent a normal operating environment, another might be the cost observed 

40 when natural gas prices are high, a third may be the price function when there is a heat wave in the area, and so on. 
One can use as many scenarios as needed so that the uncertainty in f is approximated reasonably. 

Given a set of observations, Vf v /=1 .... Jf v of the prices of utility /at time t under certain conditions s, we find the 
pricing profile that best -fits this utility using a least square approach. Here is a pseudocode of this process. 

45 . Initialization. Get a set of prices, Wf t , offered by utility /at time / under scenario s. Let )jf t be the amount of electricity 
generated by the electric utility for these price offers. 

Least Squares. The approximating function must be quadratic, convex, positive, and increasing. To find the function 
that best-fits our data while satisfying the previous constraints, we solve the following quadratic program: 



min £ t+*l:ib?. e+c / t -xi$ 2 

subject to a,* c *0. U> 

55 , b ' Z *°' 



The previous optimization model can be solved efficiently using any quadratic programming approach. See, for 
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example, R. W. Cottle, J. S. Pang, and R. e. Stone, The Linear Complementary Problem, Academic Press, 1992, and 
K. G. Murty, Linear Complementarity, Linear and Nonlinear Programming, Helderman Verlag, 1988. 

Figure 2 shows a flow chart of the previous process applied to multiple utilities, 1 n. multiple scenarios, 1 ,...,S, 

and multiple periods, 1 T. The process starts in function block 201 in which the input data is passed to this procedure. 

The input data are the historical prices, Xh s t , >1 J* , and the corresponding amounts of power generated, jri* . 

These data are usually available from the power pool and assumed to be stored on the hard disk of our computer.' ln 
many cases, such data are available through the World Wide Web (WWW). At this stage, the counter / is set to one. 
The process then moves to function block 202 in which the counter s is set to one. Function block 203 sets the value 
of / to one. The process th en proceeds to function block 204 in which the quadratic program of equation ( 1 ) is constructed 
and solved. We do not discuss the solution approach here. Most optimization packages provide specialized routines 
for handling such problems. We suggest using the Optimization Subroutine Library of I BM which provides a specialized 
function, EKKQSLV, for this operation. As a result, we obtain the values of af t , bf f and cf Jt which are used to approx- 
imate the price function, i.e., production plus profit margin, for any generation requirement. In function block 205, the 
value of t is incremented by one. Decision block 206 compares t with T. If t<T, the execution is directed to function 
block 204; otherwise, the execution proceeds to function block 207. In function block 207 the value of the scenario 
index s is increased by one. Decision block 208 checks if s is greater than S. If the answer is no, then the execution 
is directed back to decision block 203. Otherwise, the execution continues to function block 209. In function block 209, 
the value of / is incremented by one. Decision block 210 directs the flow to function block 202 if / is less than or equal 
to n, otherwise, the process terminates. As a result, we obtain a set of quadratic functions that can be used to approx- 
imate the price functions under different scenarios and during different time periods. 

The Explicit Case 

The second case we discuss is that in which the price function is known explicitly. That is, for each generating unit 
owned by utility /, we know the operating range, the cost function within that range, the start-up cost, and any other 
characteristics relevant to a generator. We may also know the status of each generator. That is, we know if a generator 
is scheduled for maintenance and if it is unavailable due to mechanical problems. Furthermore, one may have historical 
information regarding the probability of failure of each unit. All of the previous information implicitly represents the cost 
structure of a utility That is, we do not have an explicit representation of f* but we can compute the value of f* given 
the load and conditions of the generators. Note that in the current electric-power industry in the USA. the previous data 
are known to the general public. It is available on the World Wide Web and can be obtained in the form of a database 
stored on a floppy disk. However, due to competition, this information will not be available in the near future in which 
case we need to rely on the implicit representation mentioned above. 

To build a set of functions, 1 t that approximate the cost structure of a utility under different circumstances, we sort 
the generators in the increasing order of their average generating cost. This is done by adding the start-up cost of a 
unit to its operating cost at the maximum operating level and dividing by its maximum operating level. To clarify the 
previous approximation, the cost function of Figure 3A is used. To approximate this function, the start-up cost is added 
to the value of the function at the maximum operating level, G. Then, the average cost of this unit is computed by 
dividing the new cost at G by G. As a result, a generator has a linear cost function that passes through the origin and 
its operating range is anywhere between 0 and G (see Figure 3B). The generating units are sorted in the increasing 
order of their average costs which results in a convex piecewise linear function for operating the system of a utility 
Figure 4 shows an example of such a function. The first, i.e., cheapest, generator is used to cover the load between 
0 and its maximum capacity, G,. The second most economic generator is used to meet the electric load between the 
additional load up to the capacity of the second generator. That is, it covers the load between G, and G^+G^. The 
process is repeated until all the generating units of a utility are incorporated into the system. Note that we can build 
different functions to reflect the different operating conditions of a utility. For example, if some generators are not 
available, then we have a different function. Another example is the case in which a generator uses gas as its primary 
fuel. In this case, the cost function can be affected by gas prices in the market. Hence, we can create different functions 
(and different approximations) to model different market conditions. 

Note that the resulting function is an approximation of the cost function of a utility and not its price function. To 
obtain the price function from the cost function we need to estimate the profit margin of a utility. The goal of any 
deregulated market is to force the competing entities to reduce their cost. The minimum cost at which a utility can sell 
its production, without losing money, is the production cost of its electricity. However, this pricing scheme does not 
generate any profit to the utility To avoid the problem of zero profit, each competitor usually sets a profit margin that 
is acceptable for him or her. This profit margin may be different for the different periods, t, of the planning horizon and 
under the various scenarios. For example, a profit margin of 1% may be acceptable in periods in which the market 
demand is low while a higher margin might be desirable if the demand is high. The margin can also change with the 
scenarios s. For example, if fuel prices are high, one may not be able to enforce high profit margins. The price function 
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15 



is then obtained by adding a linear term that takes into consideration the profit margin resulting in ff t . 

To reflect the importance of each scenario and to make the model more complete, a probability is associated with 
each scenario s. That is, a probability, p t s , is associated with each price function that reflects the likelihood of its oc- 
currence in the future. The value of pf can be approximated using the frequency of each of the functions P. ft 
from the past. Of course, the user of the tool may adjust these probabilities using his or her intuition if there is a need 
to do so. 

3. Using the weather forecast, create a load forecast for each utility of interest. Usually, these tools require 
historical data for the load on the system of a utility and the corresponding weather conditions, i.e., maximum and 
minimum temperature during the day. Table 2 shows an example of this data. 

Table 2: 



Load vs. Weather Data 


Date 


Max 


Min 


1 


2 


3 




5/25/95 


65 


43 


2101 


2090 


1045 




5/26/95 


68 


44 


2173 


2120 


1932 




5/27/95 


74 


53 


2318 


2292 


2056 





20 



25 



35 



40 



45 



The first column is the date, second column is the maximum temperature over the twenty-four hours of that day, and 
the third column is the minimum temperature during that day. The next twenty-four columns (only three are shown) 
provide the electric load in MWH at each hour of the day. By passing this data to a forecasting tool, it can create a 
model that correlates the maximum and minimum temperatures to the load. These tools are also capable of incorpo- 
rating the effect of time of day and day of year into the forecast. The historical data of any utility can be obtained in an 
electronic form from the power pool. In the future, the electric load of a delivery point, i.e., geographic region, will be 
posted on the World Wide Web. We do not discuss the forecasting process in this document. See, for example, R. H. 
Shumway, Applied Statistical Time Series Analysis, Prentice-Hall (1988), and A. C. Harvey, Forecasting, Structural 
Time Series Models and the Kalman Filter, Cambridge (1991 ). Note that the forecast is given usually for each hour of 
the coming week. The forecasting tool also provides a probabilistic distribution, i.e., mean and standard deviation, for 
the load at each hour. That is, using the weather forecast, one can obtain for each utility a probabilistic distribution for 
its load at each hour of the planning horizon. If the forecasting tool provides a continuous distribution for the load at 
each time period, the continuous distribution can be approximated by a discrete one. The number of discrete realizations 
at time t is assumed to be K t Each demand realization, df t , i=1 has a probability of 7t t k that reflects the likelihood 
of its occurrence. 

The weather forecast is usually available on the World Wide Web and can be imported by our system easily to 
pass it on to the forecasting system. 

4. Transmission f , , in the different power lines of the electrical system. Our notation permits changing the 
transmission cost from one period to another. It also allows using multiple cost scenarios depending on the status of 
the electric system. We assume that for each transmission line (/,y) we know the maximum transmission capacity MWH 
which we denote by u* f . If a power line is not available during a certain period, we can set its capacity, f , to 
zero. Note that we define a power line by its ends or nodes / and / To simplify the task of representing the electric 
network and its delivery and generating points, we assume that there are n nodes in the network. Each node, /, is 
assumed to have a generating utility, i.e., a set of generators, associated with it and an external demand. Clearly the 
upper capacity of the generating utility, Gf f , at node / can be set to zero to model a demand center. For a generating 
node, the external load, o*, , can be set to zero. For transmittal or transshipment nodes, we set both generation and 
demand to zero. 

The transmission cost and capacity in the United States are posted on the World Wide Web as required by the 
Federal Energy Regulatory Commission. We assume that our system can access the transmission costs and capacities 
for the different regions of interest. 



Solution Approach 



55 



The computer implemented tool described herein uses the probabilistic distributions of the electric load and of the 
cost function of each utility to build probabilistic distributions for the spot-market prices and for the electric trades in 
the system. Due to the use of discrete random variables, i.e., finite S,and K h one can consider all combinations of s 
and k at each hour. The problem of cost minimization while satisfying demand and profit margins is then solved for 
each one of these combinations to produce a spot-market price and to compute the trades. The results can be bucketed 
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to form a distribution of the spot prices and trades. Given that the number of combinations at each time period is S^K^ 
it may be hard to perform all of the previous calculations in a reasonable time. Instead, one can sample from the electric- 
load and cost distributions to create a sequence of deterministic models. The sampling process (i.e., simulation) is an 
approximation of the exhaustive one described above. The sampling process is repeated until a reasonable approxi- 
5 mation is reached of the distributions of trades and spot prices or until a maximum number, L, of samples is chosen. 

Consider the market equilibrium at time / under a sample of price functions f p h<\,...,n, transmission costs c t ^ /=1,..., 
n, /=1 ,...,n, and loads d h r=1 ,...,n. We drop the indices /, s and Ho simplify notations. Due to the competition resulting 
from the bidding process, we expect the spot price to be the minimum cost at which the demand can be met while all 
suppliers satisfy their profit margins. In other words, one needs to minimize the cost of supplying power, 



n n n 

15 

The decision variable x,- represents the amount of power produced by supplier / while y t j represents the amount of 
power flowing from / to / Note that if two nodes are not connected directly with a transmission line, then the capacity 
of this line is set to zero. The minimization of (2) must be solved subject to Kirchhoff-law which requires that the total 
power flowing into a node must be equal to total power flowing out of this node. Mathematically, this condition is written 
20 as 



n n 

25 

where d- t is assumed to be zero for nodes that do not have electric load and maximum generation, G;, is zero for nodes 
that do not have any generation. To satisfy the transmission capacity, we require that O^y^Ujj, hi ,...,/?, /=1 ,...,n. Also, 
given the maximum generating capacity of utility /, we must satisfy the constraint 0<x<G h £=1,...,n. One can enforce 

30 other constraints such as the maximum capacity of a nodei 

Due to convexity, the minimization of (2) can be solved easily. See, for example, M. S. Bazaraa and C. M. Shetty, 
Nonlinear Programming: Theory and Algorithms, John Wiley & Sons, Inc., 1979. If the price functions, f h are quadratic, 
one can use a speedy quadratic programming approach as described by R. W. Cottle et at. and K. G. Murty, ibid. Note 
that by solving (2) one obtains the amount of power to be generated and transported from each supplier to each demand 

35 point. One also obtains the dual variables or the Lagrange multipliers associated with the constraints of equation (3). 
The dual variable associated with delivery point / represents the marginal cost of power at point /. That is, it is an 
estimate of the spot-market price at / under the given demand and price functions. 

We repeat the previous process by sampling a new set of price functions and electric loads then resolving the 
model. The sampling process continues until changes in the distribution of spot prices and the decision variables x,-j 

40 are within a pre-specified threshold. Another possibility is to repeat the process for a fixed number, L, of samples. 

The flow chart of Figure 5 describes the sampling process. It starts in function block 51 in which the time index, t, 
is set to 1. The process then proceeds to function block 52 in which the sample counter, I, is set to one. Note that 
whenever we select a new sample, i.e., price functions, transmission costs, and transmission capacities, the value of 
€ is incremented by one (function block 56). In decision block 53, we check if we have reached the maximum number 

45 of samples, L If the answer is yes, i.e., €=L, then the process proceeds to decision block 56 in which we move to a 
new time period in the planning horizon. If €</_, then the process moves to function block 54 in which a new set of price 
functions (one for each utility) and electric loads at the different nodes of the network is sampled. This step, function 
block 54, calls the sampling procedure of Figure 6 which returns a set of price functions, electric loads, and transmission 
capacities and costs. The data are then used to construct the problem of minimizing (2) subject to the constraints of 

so equation (3) and other system constraints. We do not discuss the solution process here. Any commercial package that 
supports quadratic programming could be used. As an output for the optimization, we obtain the amount of power 
flowing in the different lines of the network, f , and the marginal costs, n*j t , corresponding to equation (3) at the 
different nodes of the system. The process is tnen directed to function block 53 in which a new sample is generated. 
After reaching the maximum number of samples in decision block 53, the process moves to function block 56 in which 

55 the time counter, t, is incremented by 1 . In decision block 57, if f<7, then the process goes back to function block 52; 
otherwise, the process terminates. At the end of the process of Figure 5, we have for each time period a set of L values 
for the spot-market prices at nodes hi ,...,/! We also have a set of possible power flows between the different nodes 
in the transmission system. 
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Note that function block 54 of Figure 5 calls the process of Figure 6 which performs the sampling. The sampling 
process is performed using a uniform random number generator. Uniform random number generators are available on 
most computers and under any programming language. For example, the C programming language provides the func- 
tion rand as part of its standard library header. See B. W. Kernighan and D. M. Ritchie, The C Programming Language, 

s Prentice Hall, 1988. The probability of the different scenarios, pj, is used in conjunction with the electric-load proba- 
bilities, n t k f to come up with a sample. The process of Figure 6 receives from the calling process (function block 54 of 
Figure 5) the value of f. In function block 601, a random number x between 0 and 1 is sampled from the uniform 
distribution. Function block 602 defines a variable, y, and sets its value to zero. The variable y is used to compute the 
cumulative joint distribution at any given /rand s. The value of the scenario counter, s, is set to 1. Function block 603 

w sets the value of /c to 1 . In decision block 604, the value of the cumulative probability, y, is compared to x. If x is less 
than or equal to y, then the functions ff t and loads tf t are returned to the calling routine (function block 54) to be used 
in the optimization problem. Otherwise, the value of k is incremented by 1 in function block 605. Decision block 606 
checks if fchas reached its maximum limit of K t If k<Kf, then the value of the cumulative probability, y, is increased by 
p^n*" 1 and the execution is sent back to decision block 604. Otherwise, the value of s is increased by 1 in function 

15 block 608. Function block 609 updates the value of y by adding p^ 1 n** 1 to it. The execution then loops back to function 
block 603. Function block 610 is reached from decision block 604 when the value of x is less than the cumulative 
distribution at the point (k,s). 

System's Output 

20 

At termination of the simulation runs, we have the spot prices or marginal costs, X e jt at each point /, ^ t ... t n, and 
at each time period of the planning horizon. Here, the superscript £ indicates the sample index. We also obtain the 
power flow, / ^ , in the different transmission lines of the power network. Note that if / f is greater than zero, then 
node / is transmitting power to node j under sample t of our simulation. Of course, if y^ t >0 t then / =0. . To get 
25 the probabilistic distribution of the spot-market prices, X e : f , we perform the following steps'. 

Create ranges for X* jt . For example, we may want to have intervals of $5/MWH width. That is, we have the intervals 
(0-5], (5-10], (10-15] 

Count the number of X € . that falls within each range. We normalize these numbers to obtain a probability for each 
30 range. 

The resulting histogram is an approximation of the probability distribution of the marginal cost at delivery point / at 
time t. 

The previous process is called bucketing and is described in the flow chart of Figure 7. The process requires a 
35 number of buckets (i.e., ranges) which is defined by the user and the lower limit, u p of each range. We can assume 
that there is an extra range (after the last one defined by the user) that covers the range up to infinity. Function block 

71 sets the number of observations (so far) in each range or bucket to 0. The process then moves to function block 

72 in which the value of X*. f is compared against the values of uj to find the index / for which ufi X € jt <u^ +v The value of 
the counter m y is incremented by 1 in function block 73 to account for sample / Function block 74 increases the value 

40 of € by 1; i.e., we are ready to process a new sample if any. Decision block 75 checks if € is less than L If this is the 
case, then we have more samples to study and the process moves back to function block 72. Otherwise, we have 
completed classifying X € ft into the appropriate buckets. Finally, function block divides the frequency, m y by L to obtain 
the relative frequency for each range or the probability. 

We can follow similar steps to create probabilistic distributions for the amounts of power traded, f , across the 

45 different power lines, (/,/). The resulting distributions can be helpful to the decision maker in deciding the amount of 
power to be generated and its price. For a power broker, these distributions can be used to sign buying and selling 
contracts with different suppliers and consumers. 

50 Claims 

1. A computer process for use in an electricity network to manage power supply by forecasting spot-market prices 
of electric power and trading transactions at different delivery points, said method comprising the steps of: 

55 inputting maximum generating capacities of each utility competing within the market, price functions at which 

a utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost and 
transmission capacity in different power lines of the electrical system; 

using probabilistic distributions of electric load and the price function of each utility and transmission cost and 
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transmission capacity between different utilities to build probabilistic distributions for the spot-market prices 
and for electric trades in the system; and 

outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a planning horizon 
and for a plurality of delivery points to assist in said supply management. 

. The computer process of claim 1 , further comprising the step of defining different scenarios wherein the maximum 
generating capacity of each utility changes under said different scenarios. 

. The computer process of claim 1 or 2, further comprising the step of using the input weather forecast to create a 
probabilistic distribution of load forecasts for each utility of interest. 

. The computer process of any preceding claim, wherein the input price functions are given implicitly 

. The computer process of claim 4, wherein a price function of a utility is approximated by a quadratic function which 
is convex, positive and increasing. 

. The computer process of any of claims 1 to 3 wherein the input price functions are given explicitly as a cost function 
of the generating units. 

. The computer process of any preceding claim, wherein the probabilistic distributions for the spot-market prices 
and for electric trades in the system are created by sampling from electric load, price functions and cost distributions 
to create a sequence of deterministic models. 

. The computer process of claim 7, further comprising the step of minimizing the cost of operating an electric utility 
system under each sampled set of price functions, transmission costs and electric loads. 

. The computer process of claim 7 or 8 wherein the probabilistic distributions for spot prices and trades are deter- 
mined by a bucketing process comprising the steps of: 

creating ranges of spot prices; 

counting the number of spot prices that fall within each range and normalizing the counted numbers to obtain 
a probability for each range: 

and generating a histogram of the probabilities of each range as an approximation of the probability distribution 
of the marginal cost at each delivery point at a given time. 

0. A system for use in an electricity network to manage power supply by forecasting spot-market prices of electric 
power and trading transactions at different delivery points, said system comprising: 

means for inputting maximum generating capacities of each utility competing within the market, price functions 
at which a utility is willing to sell its power at a given time period, a weather forecast, and a transmission cost 
and transmission capacity in different power lines of the electrical system; 

means for using probabilistic distributions of electric load and the price function of each utility and transmission 
cost and transmission capacity between different utilities to build probabilistic distributions for the spot-market 
prices and for electric trades in the system; and 

means for outputting probabilistic distributions for spot prices and trades at a plurality of time periods of a 
planning horizon and for a plurality of delivery points to assist in said supply management. 
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